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ABSTRACT 



Context. GRB 110709B is the first source for which Swift BAT triggered twice, with a time separation of ~ 10 minutes. The first 
emission (Episode 1) goes from 40 s before the first trigger up to 60 s after it. The second emission (Episode 2) goes from 35 s before 
the second trigger to 100 s after it. [...] Within the Induced Gravitational Collapse (IGC) model, we assume the progenitor to be a 
close binary system composed of a core of an evolved star and a Neutron Star (NS). The evolved star explodes as a Supernova (SN) 
and ejects material that is partially accreted by the NS. We identify this process with Episode 1. The accretion process brings the NS 
over its critical mass, thus gravitationally collapsing to a BH. This process leads to the GRB emission, Episode 2. [...] 
Aims. We analyze the spectra and time variability of Episode 1 and 2 and compute the relevant parameters of the binary progenitor 
and the astrophysical parameters both in the SN and the GRB phase in the IGC paradigm. 

Methods. We perform a time-resolved spectral analysis of Episode 1 by fitting the spectrum with a blackbody (BB) plus a power-law 
(PL) spectral model. We analyze Episode 2 within the Fireshell model, identifying the Proper-GRB (P-GRB) and simulating the light 
curve and spectrum. We establish the redshift to be z = 0.75, following the phenomenological methods by Amati, by Yonetoku and 
by Grupe, and our analysis of the late X-ray afterglow [...]. 

Results. We find for Episode 1 a temperature of the BB component that evolves with time following a broken PL, with the slope 
of the PL at early and late times times a = and /? = -4 ± 2, respectively. The break occurs at t = 41.21 s. The total energy of 
Episode 1 is fiP = 1.42 x 10 53 erg. The total energy of Episode 2 is = 2.43 x 10 52 erg. We find at transparency a Lorentz factor 
T ~ 1.73 x 10 2 , laboratory radius of 6.04 x 10 13 cm, P-GRB observed temperature kT P - GRB = 12.36 keV, baryon load B = 5.7 x 10~ 3 
and P-GRB energy of E P _ GRB = 3.44 x 10 50 erg. [...] 

Conclusions. We interpret GRB 110709B as a member of the IGC sources, together with GRB 970828, GRB 090618 and GRB 
101023. The existence of the XRT data during the prompt phase of the emission of GRB 1 10709B (Episode 2) offers an unprecedented 
tool for improving the diagnostic of GRBs emission. 

Key words. Gamma-ray burst: individual: GRB 1 10709B — Black hole physics 

1. Introduction sion. Finally, a SN emission is either observed or expected in 

„ , . , , , . , r association with the GRB, ~ 10 days after the burst in the rest 

Of all the astrophysical processes currently being analyzed few frame Wg dm tQ find SQmces ^ which ^ data ^ of such a 

O *re more fundamental than the one presenting the coincidence d n tQ ^ tQ §ee ^ lete sequence . 

CO of some Gamma-Ray Bursts (GRBs) with the explosion of a 

— H Supernova (SN). For this, the Induced Gravitational Collapse The prototype for the IGC paradi gm has recently been given 

*^ (IGC) paradigm was first introduced by iRuffini et all J 2001b) in the analysis of GRB 090 618 (llzzo et al] 1201 21). following 

. and further analy zed in Ruffi ni et all d2007l).lRuffini et alJd2 008), the works of iRueda & Ruffinil d2012l) and llzzo, Rueda & Ru ffini 

^ IRueda & Ruffinil (120121) and llzzo. Ru eda & Ruffini ( 20121) . Re- d2012l) . In this work we follow the same line and identify four 

H cently, it has been evidenced that indeed this process can explain different episodes in GRB 110709B. Episode 1 starts 40 s be- 

. . .the coincidence between the SN and GRB emis sion, both from fore the first trigger and lasts up to 60 s after it and is well fit 

an observational and a th eoretical point of view (llzzo et alJ2Q12b by a blackbody (BB) plus a power-law (PL) spectral model . 

IPenacchioni et alj|2012l) . It coiTesponds to the trigger of the SN explosion of the com- 

In the IGC paradigm (Ruffi ni et alJl200fbl 120071) . a binary pact core and its accretion onto the NS companion. We notice 

system formed by an evolved star and a Neutron Star (NS) com- that the B B temperatur e decays with time following a broken 

panion is considered as the progenitor. power-law (Ryde 2004). Episode 2 starts 35 s before the second 

The IGC paradigm implies a well determined time sequence, trigger and lasts up to 100 s after it. It corresponds to the emis- 

In a close binary system of a massive star in the latest phases of sion of the canonical GRB emitted in the formation of a BH. 

its thermonuclear evolution and a NS companion, the massive Episode 3 starts at 800 s all the way to 10 6 s. It consists in a 

star undergoes a SN explosion. The accretion of the early-SN standard X-ray emission identified in all systems following the 

material onto the NS companion leads the NS to its critical mass IGC paradigm (Pisani et al., in preparation). Episode 4 corre- 

and consequently to its gravitational collapse to form a black sponds to the observation of the optical SN emission, observable 

hole (BH). The emission of a canonical GRB in the collapse to after T b S — (1 + z)Tsn- In the present case, there is no evidence 

the BH takes place. A young NS is born out of the SN explo- of an associated SN in the optical band. An explanation for this 
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is given bv lZauderer et al.1 (120121) . who classified GRB 1 10709B 
as dark and stated that its optical emission may have been ab- 
sorbed by the host galaxy and/or the interstellar medium. The 
ensemble of these four episodes characterize the IGC scenario. 

As an outcome, at the endpoint of the IGC scenario, a binary 
system represented by a NS (formed by the SN explosion) and a 
BH (formed after the GRB explosion) should be expected. 

As in the case of GRB 101023, we do not know the cos- 
mological redshift of GRB 110709B due to the lack of optical 
data. Therefore, we in fer it from ph enomenological methods: 
1) The Amati re lation (|Amati 20 0g), 2) the Yoneto ku relation 
dYonetokul2004llbl), 3) the workoflGrupe et alj d2007b and 4) the 
work bv IPenacchio ni et al.l d2012hJRuffinil (1201 2l) and by Pisani 
et al. (in preparation), which describe a scaling of the late X- 
ray emission of GRB 090618. In the case of GRB 111228, 
which we are currently analyzing, we find a striking coincidence 
between the values of the cosmological redshift determined by 
these methods for GRB 1 1 0709B . 

In section [2] we report the observations of the two compo- 
nents of GRB 110709B by the different instruments, in space 
and on the ground. In section[3]we reduce the Swift data and per- 
form a detailed spectral analysis of both Episode 1 and Episode 
2. In section |4] we infer the redshift of the source using the four 
phenomenological methods mentioned above. In section [5] we 
determine the radius of the emitting region from the knowledge 
of the redshift and the BB flux of the first episode. In section [6] 
we give a brief description of the Fireshell model and we per- 
form a deeper analysis of Episode 2 within this model, repro- 
ducing the light curve and the spectrum by a numerical simu- 
lation. In section [7] we calculate the parameters of the binary 
progenitor leading to the IGC of the NS to a BH by the SN ex- 
plosion. Details on the accretion rate onto the NS, total accreted 
mass, SN ejecta density, NS mass, and binary orbital period are 
obtained for selected values of the SN progenitor mass. In sec- 
tion [§1 we comment on th e radio emission detected by EVLA 
dZauderer & Bergerll2012l) . In section we present the conclu- 
sions. 

2. Observations of GRB 110709B 

GRB 1 107 09B has been detected by Suzaku dOhmorietal.l201ll) 
and Swift (Cummings et al.l 1201 1[) satellites, and by ground- 
based telescop es like GROND dUpdike et aljfeoi lh and Gemini 
(Berg eril201ll) . 

The Burst Alert Telescope (BAT) on board Swift triggered a 
first time at 21:32:39 UT (trigger N°= 456967). The location of 
this event is RA= 164.6552, DEC= -23.4550. The light curve is 
composed of multiple peaks, with the whole emission extending 
up to 60 s after the trigger (see Fig. [T}. What is most inter- 
esting is that there was another trigger at 21:43:25 UT (trigger 
N°=456969), ~ 11 minutes after the first trigger. The on-board 
calculated location is RA=164.647, DEC=-23.464. This time 
Swift did not need to slew, because it was already pointing to 
that position. This second emission shows a bump that begins 
100 s before the second trigger and lasts around 50 s, followed 
by several overlapping peaks with a total duration of about 40 s, 
and another isolated peak of 10 s of duration, 200 s after the sec- 
ond trigger. Fig. [T]shows the complete BAT light curve and Fig. 
|2] shows the light curve taken by the X-Ray Telescope (XRT) in 
the 0.3-10 keV band. 

There have not been detections in the optical band by Swift 
UVOT, which started to observe 70 s after the first BAT trigger 
dHollandll20i"ll). The observat ions with GROND at La Silla Ob- 
servatorv dUpdike et al.ll20TTl) simultaneously in the g' r' i 'z' 




Fig. 1. BAT light curve of GRB 110709B, including both triggers. 
Here we can appreciate the time separation (about 10 minutes) between 
the first and the second trigger. The light curve is in the (15-150 keV) 
energy band. The time is relative to the first trigtime, of 331939966 
s (in MET seconds). The second trigger was at 331940612 s in MET 
seconds. 



XRT data of GRB 110709B 
Cyan: WT setting; Blue: WT; Red: PC 
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Fig. 2. Count light curve of GRB 1 10709B obtained from the Swift- 
XRT detector, in the (0.3-10 keV) energy band. The time is relative to 
the first trigtime, of 331939966 s (in MET seconds). 



JHK, reveal tw o point sources within the 5" .3 XRT error cir- 
cle reported bv ICummings et alj (1201 lh . They suggest that one 
of them could be an afterglow candidate for GRB 1 10709B, al- 
though it is very faint. 

It has been suggested by Zaud erer et al.l (1201 2l) that this 
source is an "optically dark" GRB. The possible reasons for this 
are: 1) dust obscuration, 2) intrinsically dim event, and/or 3) 
high redshift (optical emission suppressed by Lya absorption at 
A b s < 1216A(1 + z)). However, they rule out the possibility 
of a high redshift event due to the association with an optically 
detected host galaxy. Furthermore, they have inferred the opti- 
cal brightness of the after glow according to the standard after - 
glow synchrotron model (Granot&Sari 20021 ISari et al .1 119991) . 
and from the non-detection in the optical-NIR wavelengths they 
find a very large rest frame extinction for GRB 1 10709B . This 
can explain the lack of detections in the optical band. 

There have bee n detections in the radio band on several oc- 
casions by EVLA (Zauderer & Berger 2012), revealing a single 
unresolved radio source within the XRT error circle, which re- 
brightened by a factor of 1.6 between 2.1 and 7 days after the 
burst. The location of the source is RA=10:58:37.114, DEC=- 
23:27:16.760. 



Article number, page 2 of 1 1 2 1 



A.V. Penacchioni et al.: GRB 110709B 



Unfolded Spectrum 



GRB 1 1 0709B light curve (1 st trigger) i 



mi 



1 1 a 



Fig. 3. Count light curve of Episode 1 of GRB 1 10709B obtained from 
the 5w(//-BAT detector, in the (15-150 keV) energy band. The time is 
relative to the first trigtime, of 331939966 s (in MET seconds). 



3. Data analysis 

In the following we refer to the emission that goes from 40 s 
before the first BAT trigger to 60 s after it as Episode 1 (see Fig. 
0. We call the emission going from 35 s before the second BAT 
trigger to 100 s after it as Episode 2. We make use of Swift BAT 
data to perform the spectral analysis with XSPEC. 



3.1. Episode 1 

We perform a time-integrated analysis to the whole Episode 

1, using five di fferent spectral models, namely BB, Band 
dBandetal.il 19931) . BB+PL, PL and CutoffPL. The results of the 
fits are shown in Tableland in Fig. |4] The Band function is not 
well constrained, so we have excluded it in the following anal- 
ysis. A statistical test shows that the best models are BB+PL 
(X 1 = 56.65) and CutoffPL Of 2 = 54.45). Since the difference 
in the x 2 between these two models is 2.2, the two models are 
statistically equivalent. So we discriminate between these two 
models based on the physical grounds expected from the IGC 
scenario. In this scenario, we expect a thermal emission from 
the expansion of the outer layers of the compact core SN pro- 
genitor. Thus, we have chosen the BB+PL model. We obtain a 
BB temperature kT = (22 ± 5) keV, a PL index y = 1.4 + 0.1 
and a x 2 = 56.65 (54 DOF). The flux of the BB component 
is ~ 12% of the total flux. The total energy of Episode 1 is 
Ejsa ~ 1-42 x 10 53 erg. The results of the fit are shown in Table 
Q] Then we perform a time-resolved spectral analysis with a bin- 
ning of 5 s fitting the same model and find that the temperature 
of the BB comp onent follows a broken power-law, as mentioned 
in Ryde (2004), from 5 s before the trigger to 55 s after it (see 
Fig. 0. The broken power-law is indeed a constant function 
plus a simple power-law function. This is t he same behav ior as 
for the previo usly analyzed GRB 0906 18 (llzzo et al.ll2012l) and 
GRB 101023 dPenacchioni et alJl2012l) . However, we notice that 
the temperatures for this GRB are lower. Nevertheless, the si- 
multaneous presence of a BB and PL component is necessary in 
order to obtain an acceptable fit of the data (see Fig. |4). 

3.2. Episode 2 

We also performed a time-integrated spectral analysis of Episode 

2, whose light curve is shown in Fig. Q This episode starts 35 
s before the second trigger and last 135 s, until 100 s after the 




Fig. 4. Fit of Episode 1 with a BB+PL model. The parameters of 
the fit are: kT = (22 ± 5) keV, BB Amp= 0.2 ± 0.1, y = 1.4 ± 0.1, PL 
Amp= 2.2 ± 0.8, Red-^ 2 = 1.049 (54 DOF). 




Fig. 5. Evolution of the kT component of the BB+PL model during 
Episode 1. The first data point corresponds to 5 s before the first BAT 
trigger. The vertical line corresponds to the trigger time. The time is 
in the observer frame. The temperature evolves in time following a 
broken power-law fit. There is a break at t = 41.21 s. The indices of 
the PL are a = (consistent with a constant function) and /3 = -4 ± 2, 
respectively. The presence of the BB, although smaller than in previous 
cases, is essential to have an acceptable fit. 



second trigger. We tried to fit the spectrum with the following 
spectral models: BB, PL, BB+PL, cutoffPL and Band (see Ta- 
ble We can easily discard the BB and Band models because 
in one case the Red x 2 is too high and in the other case there 
is an unconstrained parameter. As the PL and the BB+PL are 
nested models, we performed a statistical test to see which one 
is the best. We obtained a probability Prob=0.001 that the sim- 
pler model is better, so the BB+PL dominates over the PL. Then 
we have to compare between the BB+PL and the CutoffPL mod- 
els. As they are not nested, we cannot apply the same test. So 
we chose the model that gives the lowest^ 2 . We concluded that 
the model that best fits Episode 2 is the cutoffPL model. 

It is clear from the analogies with GRB 090618 and GRB 
101023 that Episode 2 has all the characteristics of a canonical 
GRB. A difference between GRB 1 10790B and the already ana- 
lyzed ones is that the separation between Episode 1 and Episode 
2, ~ 10 min, is much bigger than previously, ~ 50 s. This re- 
markable time separation between the two episodes is an ad- 
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Table 1. Fit results of Episode 1 with five spectral models: BB, Band, BB+PL, PL and CutoffPL. The flux is in the energy band (15- 150) keV, in 
units of erg/cm 2 /s. 



DD 


RnnH 
DallLl 


1)1) I 1 , 


pi 


PntnfTPI 


kT = 18.9 ±0.2 


a = -1.2 + 0.1 


kT = 22 ± 5 


y = 1.37 + 0.02 


y = 1.1 +0.1 


K BB = 0.95 + 0.01 


/3 =unconstr. 


Kflfl = 0.2 + 0.1 


K P0 = 2.0 + 0.2 


E = 196 + 68 




E = 296 + 255 


y = 1.4 + 0.1 




K = 0.8 + 0.2 






K P0 = 2.2 + 0.8 






Redx z = 7.3 


Red x z = 1.031 


Redy 2 = 1.049 


Redx 2 = 1.14 


Red* 2 = 0.99 


56DOF 


54 DOF 


54 DOF 


56 DOF 


55 DOF 


Flux= 7.52 x 10~ 8 


Flux= 8.99 x 10~ 8 


Flux= 8.96 x 10~ 8 


Flux= 9.08 x 10~ 8 


Flux= 8.93 x 10~ 8 



Unfolded Spectrum 




Fig. 6. Fit of Episode 2 with a Cutoffpl model. The photon index is 
y = 1.0 + 0.1, the normalization constant is 0.5 ±0.1, the cutoff energy 
is E = 132 ±31 and the reduced chi squared of the fit is Red-^ 2 = 0.77 
(55 DOF). 




50 200 250 



Fig. 7. Count light curve of Episode 2 of GRB 1 1 0709B obtained from 
the 5wi/f-BAT detector, in the (15-150 keV) energy band. The time is 
relative to the second trigtime, of 331940612 s (in MET seconds). 



ditional new fact to propose a different astrophysical origin of 
these two components. 

We turn now to the crucial analysis of the determination of 
the cosmological redshift of Episode 2. 



4.1. N H column density 

We first tried to ge t an upper limit for z following the work of 
Grupe et al. (2007). They consider a relation between the ab- 
sorption column density in excess of the galactic column density, 
given by ANh - Nhjh - Nn, ga i and the redshift z, through the 
equation 



log(l +z) < 1.3 - 0.5 [log(\ + AN H )]. 



(1) 



We calculated Nn,gai from the radio map of the galaxy in the 
Lab Survey websit^B by entering the coordinates of the GRB 
(RA=164.64, DEC=-23.46). We obtained N H , ga i = 10.5 X 10 20 

-2 



cm 



To obtain the value of Nhju we took the XRT data from 
3 2Q0fl»s to 10 6 s after the first BAT trigger and fitted the model 
phabs*po using the program XSPEC. The XRT data were re- 
duced by the xrtpipeline software, version 0.10.4, which is part 
of the HEASOFT package, version 6.12. We use the standard 
response matrix swxpc0tol2s6_20010101v013.rmf for the PC 
mode data. The model phabs represents the photoelectric ab- 
sorption 



M(E) = e - n '" T(E \ 



(2) 



where n# is the equivalent hydrogen column density (in units 
of 10 22 ctrr 2 ) and cr(E) is the photoelectric cross section, not 
including Thompson scattering. We obtained a value of Nhju - 
71.76 x 10 20 cirT 2 . Using these values in ([1) we obtained an 
upper limit for the redshift of z < 1 .35. 



4.2. Amati Relation 

We also tried to determine th e redshift of Episode 2 through the 
Amati relation (lAmatil 12006). that relates the isotropic energy 
Ei S0 of the GR B to the peak Ener gy in the rest frame E p j of the 
vF v spectrum (Ama ti et al.ll200 9). The analytical expression of 

Ejso is 



(1+z)' 



<bol, 



(3) 



where d\ is the luminosity distance, z is the redshift and Sboi is 
the bolometric fluence, related to the observed fluence in a given 
detection band (£ m ,„, E„ mx ) by 



4. Cosmological redshift determination 

We used four phenomenological methods to constrain the red- 
shift of the source, based on different relations, detailed below. 



r 10 4 /(l+z) , 

Ji/(i + z) E <t>(E)dE 



» bol 



* obs ~ 



JJT E4>(E)dE 



(4) 



http://www.astro.uni-bonn.de/~webaiub/english/tools_labsurvey.php 
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Table 2. Spectral fit of the whole Episode 2 with different models: BB, PL, BB+PL, cutoffPL and Band. The flux corresponds to the (15 -150) 
keV energy range. The model that best fits the data is the cutoffPL. We can easily see that the BB and Band models are not good. As the PL and 
BB+PL are nested models, we performed a statistical test to see which was the best one. We obtained a probability Prob=0.001, indicating that 
the BB+PL dominates over the PL model. Then, between BB+PL and CutoffPL (that are not nested, then we cannot apply the same test), we took 
the CutoffPL as the best model because it gives the lowest^- 2 . Note that the models with which we fit the data are those defined in the XSPEC's 
Manual: http://heasarc.nasa.gov/xanadu/xspec/xspecl 1/manual/manual.html 



BB 


PL 


BB+PL 


CUTOFFPL 


BAND 


kT [keV]= 17.5 ± 0.2 
Kbb — 0.661 + 0.009 


y= 1.46 + 0.02 
K P0 = 2.1 +0.2 


kT [keV]=20 + 3 
Kbb — 0.16 + 0.06 
y = 1.5+0.1 
K P0 = 2.3 + 0.8 


y= 1.0 + 0.1 
E = 132 + 31 
£ = 0.5 + 0.1 


a = -1.0 + 0.1 

y8 = unc 
E = 142 + 42 
K = 0.0048 + 0.0008 


Red^ 2 =7.16 
DOF= 56 
Flux= 5.2 x 10~ 8 


Red^- 2 = 1.109 
DOF= 56 
Flux= 6.52 x 10~ 8 


Red x 2 = 0.78 
DOF=54 
Flux =6.35 x 10~ 8 


Red x 1 = 0.77 
DOF= 55 
Flux= 2.43 x 10~ 8 


Red x 2 = 0.79 
DOF= 54 
Flux=6.36x 10~ 8 



Here, cf> is the spectral model considered for th e spectral data 
fit; in this case a Band model (B and et a"flll993h . composed of 
two smoothly connected power-laws. E p j is related to the peak 
energy E p in the observer frame by 



EJl+z). 



(5) 



The peak energy is the energy at the peak of the vF v spectrum. 
It can be written as 

E p = Eq(2 + a), 

where Eo is the energy at which the two power-laws intersect 
and a is the slope of the low-energy power-law, according to the 
Band model. 

We calculate the luminosity distance di, as given by the stan- 
dard cosmological model 



«o Jo 



dx 



VQ,„(i + -*) 3 + a A 



(6) 



where the Hubble constant is Ho = 70kms~' Mpc _I ,Q,„ = 0.27, 
Qa = 0.73 and c is the speed of light. 

Following th e same procedure as described in 

dPenacchioni et alJ l2012h . we calculated Ej S0 and Epj for 
different values of z, from 0.1 to 3, at steps of 0.1. Fig. [HJshows 
that the relation is satisfied for values of z > 0.4. This puts a 
lower limit to the estimation of the redshift. 



4.3. Yonetoku Relation 

We finally obtaine d a range of po ssible redshifts by using the 
Yonetoku relation (Yonetoku 2004). This relation, also known 



as the E p - Luminosity relation (E p - L) , connects the observed 
isotropic luminosity L in units of 10 52 erg s _1 with the peak en- 
ergy E p (l + z) in the rest frame of the GRB. It is valid for val- 
ues of E p between 50 and 2000 keV, and a luminosity range of 



10 50 - 10 54 ergs-'. 

The best fit function for the E„ 



L relation is 



10 52 erg s 



(2.34!? ^) x 10 s 



E P (l+z) 



IkeV 



2.0±0.2 



(7) 



The peak luminosity and the peak energy are calculated by 
integrating within a 1 s interval around the most intense peak of 
the light curve, because this is a better distance indicator than 
the burst average luminosity. However, we took a 10s interval 




Fig. 8. Amati relation (solid line) with its 1 - <r uncertainties (doted 
lines) and peak energy E P j vs. E iso for GRB 110709B, for different 
values of the redshift, from 0.1 to 3, at steps of 0.1. We can see that the 
data matches the theoretical function within 1 - cr for z > 0.4. 



around the most intense peak in order to better constrain the 
value of the parameters (i.e., to increase the number of photons 
in the spectrum and obtain an error which is smaller than the 
value of the parameters). The peak luminosity in the rest frame 
(with the proper K-correction) can be calculated as 



L = 4nd 2 L F bo i, 
where 



E bol — *abs 



(-10000/(1+2) 
Jl/(l+z) 



EN(E)dE 



JT "'" N(E)dE 



(8) 



(9) 



is the energy flux and P„i, s is the photon flux. 

Fig. |9]shows the Yonetoku relation (solid line) with its uncer- 
tainties (dotted lines), and the values of L and E p j for each value 
of z, from 0.1 to 3, at steps of 0.1. We see that the Yonetoku 
relation is satisfied within lcr for values of the redshift > 0.7, 
consistent with the results obtained with the Amati relation. 

In conclusion, if we put together the three methods, we have 
a range of possible redshifts of 0.7 < z < 1.35. 



4.4. Estimate of the redshift using the X-ray afterglow 

We already presented in dPenacchioni et al.ll2012l) a method to 
estimate the redshift of GRB 101023 by comparing its X-ray 
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Fig. 9. Yonetoku relation (solid line) with its 1 — er uncertainties (doted 
lines) and peak Luminosity vs. E pi for GRB 110709B, for different 
values of the redshift, from 0.1 to 3, at steps of 0.1. We can see that the 
data matches the theoretical function within 1 — cr for z > 0.6. 



light curve to the one of GRB 090618, of known redshift (z — 
0.54). Here we rescale the X-ray light curve of GRB 090618 as 
if it was seen at different redshifts and plot it together with GRB 
1 10709B light curve, looking for the values of z for which these 
light curves overlap at late times. We find a remarkable consis- 
tency between this method and the phenomenological methods 
already mentioned. 

In oder to compare in a common rest frame the two emissions 
from the GRBs, we apply the following operations only to GRB 
090618: 

1) determination of the starting time T slarl of the late decay 
emission, 

2) spectral analysis of this emission with an absorbed power- 
law model, 

3) extrapolation of this spectral model in a common cosmo- 
logical rest-frame energy range and, consequently, rescaling of 
GRB 090618 light curve for the different energy ranges, 

4) cosmological correction for the arrival time by taking into 
due account the different scaling due to cosmological redshift, 
and 

5) correction of the observed flux by changing the redshift of 
GRB 090618. 

A detailed description of the method will be given in a forth- 
coming publication (Pisani et al., in preparation). 

In this way we compare directly both light curves for differ- 
ent redshifts of GRB 090618. Fig. [I0]shows GRB 090618 light 
curve seen as if the source was located at different redshifts: 0.2 
(blue), 0.4 (green), 0.7 (grey), 1.0 (orange) and 2.0 (purple). The 
red light curve corresponds to GRB 110709B. We can see that 
it lies between the green and the orange ones. A more accurate 
scaling of the late X-ray afterglow suggests a redshift of z — 0.75 
for this source. 

Fig. [TQshows the superposition of GRB 1 10709B and GRB 
0906 1 8 light curves in the observer frame, as if they were loc ated 
at a redshift z = 0.75. 

There is however a second aspect which is due to the pecu- 
liarity of the turn-on 7q of the XRT detector. At the time BAT 
triggered for the second time, XRT was already pointing at the 
source and was able to detect the emission at very early times, 
making this GRB probably the first for which XRT has the ear- 
liest detection up to date. We need to shift GRB 1 10709B light 
curve in order to make the early steep decays (originating in the 
prompt in our interpretation) coincide. This is done by adding 
a time T* = +800 s to GRB 090618 light curve. The superpo- 
sition is very good. In this way we also make the early decays 
coincide. This factor is arbitrary, but we need to include it be- 



cause GRB 110709B XRT light curve presents many spikes at 
the beginning, which according to our interpretation correspond 
not to the steep decay of the X-ray light curve but to the prompt 
emission. 

In the case of GRB 1 10709B, thanks to the fact that XRT was 
already active and collecting data at the time of the second BAT 
trigger, we were able to follow the behavior of the whole GRB 
emission of Episode 2. This is a key point to our understand- 
ing of GRB 110709B, since only in very few cases XRT had a 
response during the early emission. 

In the big flare at ~ 1000 s after the first BAT trigger we 
notice a strong correlation between the emission in X-rays and 
in y-rays. We identify this emission as the prompt emission of 
Episode 2. After this prompt phase the traditional plateau phase 
is observed. After the plateau phase, there is the late decay phase 
in the X-ray light curve following a power-law behavior which 
has already been observed in other sources (i.e., GRB 101023, 
GRB 090618, GRB 111228). We study this decay in the IGC 
paradigm, and consider the possibility that it is produced by the 
early emission of the newly-born NS. It is interesting to notice 
that in GRB 110709B the typical flare in X-rays just preced- 
ing the plateau phase and following the prompt emission is not 
observed. This X-ray emission usually occurs without any asso- 
ciated y-ray emission, since the data is usually below the BAT 
threshold. In the present case, it is conceivable that the flaring 
indeed occurred during some of the gaps of ~ 4000 s in which 
there is no data due to Earth occultation. 

We can then distinguish two types of flares in the X-ray 
light curve. The first type occurs at early times, previous to the 
steep decay, and belongs to the prompt emission. This flares can 
be seen in X-rays only when XRT starts its detection at early 
enough times, e.g., when the satellite was already pointing at 
a region near the burst position and did not need much time to 
slew. The light curve in X-rays generally follows the trend of 
the light curve in y-rays. The second type of flares occurs at 
later times, just preceding the plateau phase. This flares are seen 
only in X-rays since their photon flux is much lower than the 
BAT threshold. In the ICG paradigm, we interpret this flares as 
possible indicators of the breakout of the SN. 

We are currently analyzing more sources in the catalogue by 
Mar gutti et al.l (1201 2|) to look for these three very distinct phases, 
i.e., the flares in the prompt emission, the flares in the afterglow 
and the late decay after the plateau, each of them having a dif- 
ferent physical origin within the IGC paradigm. 

5. Episode 1 : radius of the emitting region 

With the knowledge of the redshift and the parameters of the fit 
with a BB + PL model, we computed the isotropic energy of the 
whole Episode 1, = 1.42 X 10 53 erg. 

With the energy flux of the BB component 4>bb as a function 
of time from the time-resolved spectral analysis and the lumi- 
nosity distance di, we can compute the value of the radius of the 
emitter in cm (we then express it in km in Fig. [T2t through 



"BB 



(10) 



o-r 4 (i+z) 2. 

Here <Pbb is the BB flux in units of erg cm~ 2 s~ l , cr = 
5.6704 erg cm 2 s~ l K~ A is the Stefan-Boltzmann constant and di 
is the luminosity distance in cm. 

The best fit of the expanding radius is 



r(t) = at 



(ID 
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Fig. 10. Plot of GRB 090618 seen as if it was at different redshifts: 
0.2 (blue), 0.4 (green), 0.7 (grey), 1.0 (orange) and 2.0 (purple). The 
red light curve corresponds to GRB 110709B. We can see that it lies 
between the green and the orange light curves, indicating that the red- 
shift must be between 0.4 and 1.0. 
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Fig. 11. GRB 090618 light curve (blue) as if it was seen at z = 0.75 
together with GRB 1 10709B light curve (red). There is an excellent su- 
perposition in the late decay assuming a temporal shift of GRB 090618 
of T, = 800 s to the right to have the coincidence of the steep decays of 
the light curves. 

where a = (1.5 + 1.2) x 10 4 km s~ b and b = 0.32 + 0.27 (see Fig 

We associate the BB component to the expansion of the 
ejected material, while the power-law is associated (as we in- 
terpret from the IGC paradigm) with the accretion of part of this 
material onto the NS companion. 

6. Analysis of Episode 2 in the Fireshell model 

We recall that the Fireshell model (IDamour & Ruffinil [l97l 
IRuffini et al.ll2000t lRuffini]l200lt iRuffini et al.ll2010bl) is an al- 
ternati v e to the Fireball m odel, first propos ed by Cava llo & Reesl 
dl978l) . iGoodmanl (119861) and iPaczvnskvl il986|). We assume, 
within the Fireshell model, that all GRBs originate from the 
gravitational collapse of a star approachin g asymptotically the 
forma tion of a Kerr-Newmann BH ( Wilts hire. Visser & Scott I 
2009). An electric field E is created just outside the collaps- 
ing core a nd in between the expanding ou ter shells that act as a 
capacitor (Preparata, Xue & Ruffinil [T998I) . This electrical field 
grows until it reaches a critical value, E c - m 2 c 3 /he. At this 
time, vacuum polarization occurs, leading to pair creation at 
the expenses of the gravitational energy. An optically thick e ± 
plasma forms with total energy E e m in the range 10 49 - 10 54 erg. 
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Fig. 12. Radius of the emitting region as a function of time (in the 
cosmological rest frame), corresponding to Episode 1. The radius in- 
creases with time following a power-law at b , with a = (1.5 ± 1.2) x 10 4 
km/s and b = 0.32 ± 0.27. 



The e ± plasma reaches thermal equilibrium on a timescale of 
10~ 12 s {\k senov. Ruffini & Vereshch agin 2007). Being opti- 
cally thick, the plasma self-acceler ates due to its internal radi- 
ation pressure (Ruf fini et al.| [l999a b). After an early expansion 
in vacuum, the e ± -photon plasma engulfs the baryonic matter 
Mb of the outer shells and reaches thermal equilibrium with it. 
The baryonic matter is described by the dimensionless parame- 
ter B = Mbc 2 IE e tot . B must be less th an 10" 2 , otherwise there 
will not be any relativistic expansion (Ruf fini et alJuO OO). The 
optically thick fireshell composed by e ± -photon-baryon plasma 
self-accelerates to ultrarelativistic velocities, finally reaching the 
transparency condit ion. A flash of r adiation is then emitted. This 
is the P-GRB (Ruf fini et al.ll20 01a). The amount of energy radi- 
ated in the P-GRB is only a fraction of the initial energy E e tot . 
The remaining energy is stored in the kinetic energy of the op- 
tically thin baryonic and leptonic matter fireshell that expands 
ballistically and starts to slow down due to the inelastic col- 
lisions with the Circumburst Medium (CBM). This interaction 
gives rise to a multi-wave length emission, the extended after- 
glow ( Ruffini et al . 2001a). We can estimate the characteristic 
inhomogeneities of the CBM by fitting the luminosity of the X- 
ray source and imposing the fully radiative condition in the col- 
lision between the ultra relativistic baryonic shell and the clouds 
of the Interstellar Medium (ISM). The compl e te analytic so lu- 
tion has been developed in Bianco & Ruffini (2004, l2005allbl) . 
together with the analytic expression of the Surfaces of Equal 
arrival Time of the photons at the detector (EQTS). The after- 
glow presents three different regimes: a rising part, a peak and a 
decaying tail. We therefore define a "canonical GRB" light curve 
with two sharply different components: 1) the P-GRB and 2) the 
extended afterglow. What is usually called "Prompt emission" 
in the current GRB literature mixes the P-G RB with the raisin g 
part and the peak of the extended afterglow (Ru ffini et alJ 2003). 
The spectrum of the extended afterglow is initially assumed to be 
thermal in the comoving frame of the expanding shell. Recently, 
after the analysis of some highly energetic sources observed by 
Swift and Fermi satellites, this assumption of a pure comov- 
ing thermal spectrum has been relaxed and a phenomenological 
modification by a power- law of the low energy spectral slope has 
been introduced (Patric eli et al. 1 120121) . The observed non ther- 
mal spectrum shape is due to a double convolution of thousands 
of instantaneous comoving spectra, with different temperatures 
and different Doppler factors, ove r both the EQTS and the ob- 
servation time (Ru ffini et alll2004bl) . 
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Table 3. Fit of the P-GRB and the afterglow of GRB 1 10709B, Episode 
2. The P-GRB is well fit with a BB model, while the whole Episode 2 is 
best fit by a CutoffPL model. From this fit and the value of the redshift 
we are able to calculate E im and E P -grb- 



Parameter 


P-GRB 


P-GRB + Afterglow 


kT [keV] 


14+ 1 




BB Amp 


0.30 + 0.02 




7 




1.03 + 0.1 


PL Amp 




0.5 + 0.1 


Red^ 2 


1.448 (56 DOF) 


0.77 (55 DOF) 


Energy Flux 


2.413 x 10~ 8 


6.34 x 10~ 8 


(15- 150 keV) 






[erg crrT 2 s -1 ] 






Energy [erg] 


3.44 x 10 5l) 


2.43 x 10 52 



Having fixed the value of the redshift to z — 0.75, we started 
the analysis of Episode 2 within the Fireshell model. We first 
looked for the P-GRB during the first bump of Episode 2 (from 
100 to 40 s before the second trigger) by fitting the data with 
a BB + PL model. We selected several time intervals as the P- 
GRB during the first bump of Episode 2, but in some cases the 
fits were not good. In some other cases, to reproduce the ratio 
between the P-GRB energy and the total energy we needed to 
consider a baryon load B > 10~ 2 (which has no sense within 
the Fireshell model) and, in other cases, there was a discrep- 
ancy between the observed temperature and the one given by the 
simulation. Thus we concluded that this bump should belong to 
Episode 1 . The reason why we do not find a strong thermal sig- 
nature in this bump is that Episode 1 starts ~ 10 minutes before 
the beginning of the bump and the temperature of the BB com- 
ponent decreases very rapidly following a power-law in the first 
seconds of emission. Consequently, after such a long time we do 
not expect to find any signature of a BB from Episode 1 . 

We finally selected the P-GRB as the 9 s from 35 to 26 s 
before the second trigger, and the following emission from -26 
to 100 s as the afterglow. Table|3]shows the parameters of the fit. 
We calculated a P-GRB energy of E P - CRB = 3.44 x 10 50 erg and 
an isotropic energy of E/ so = 2.43 x 10 52 erg. 

We inserted these values of the energies into our numerical 
code and calculated the value of the baryon load, B = 5.7 x 10~ 3 . 
We simulated the light curve and the spectrum, obtaining, at the 
transparency point, a Laboratory Radius r tr = 6.04 x 10 13 cm, 
a gamma Lorentz factor F = 1.73 x 10 2 and a P-GRB observed 
temperature (after cosmological correction) kT = 12.36 keV. 

Figs. [T3b and[T3b show the simulation of the light curve and 
the spectrum of Episode 2, respectively. The photon index of 
the XRT and BAT spectra are in agreement with that predicted 
by the simulation. Details of this calculation will be given in a 
forthcoming letter (Penacchioni et al., in preparation). Fig. [T4l 
shows the density mask of the ISM, i.e., the density of particles 
of the interstellar clouds as a function of the distance to the cen- 
ter of the BH. This density has to be interpreted as an effective 
density because fragmentation may occ ur in the expanding shell 
dRuffini et alJl2007l:Dainotti et alJl2007h . 

7. Nature of the Progenitor 

Following the wor k s of iRueda & Ruffinil d2012[) and 
Izzo, Rueda & Ruffini (2012), we suggest for the origin of 
GRB 110709B a binary system formed by a massive evolved 
star on the verge of a SN explosion and a NS. The early-SN 
material expanding at non-relativistic velocities is then accreted 
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Fig. 13. Simulation of the BAT (a) light curve and (b) spectrum of 
Episode 2. We included XRT data in the fit of the spectrum to show that 
the slope predicted in the fireshell model is in agreement with the slope 
of the X-ray spectrum. 



n 



Fig. 14. Particle density of the ISM clouds as a function of the dis- 
tance. The mean density is 76 part/cm 3 . 

by the NS companion at times larger than forcer, when the ma- 
terial reaches the NS gravitational capture region. The emission 
observed in Episode 1 is associated to this early-SN evolution, 
identified with the thermal component, and accretion process 
onto the NS, possibly related to the non-thermal component. 
The NS reaches in a time fo.accr + Af accl - the critical mass and 
gravitationally collapses to a black hole, emitting the GRB seen 
in Episode 2. We assume t he critical mass of a no n-rotating NS 
M cri t = 2.67M as given bv lBelvedere et all (120121) . 

The amount of material that reaches the NS gravitational 
capture region 

s , (0 = ™ (12) 

per unit time is given by (see IRueda & Ruffinil d2012l) and 
llzzo. Rueda & Ruffinil d2012l) ) 

M(t) = npz)(t)v eiM (t)R 2 cap (t), (13) 
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where /? cap is measured from the NS center. 

In these expressions, p e j(?) = 3M e j(?)/(47rr^(f)) is the density 
of the ejecta, 

M ej (f) = M ej (0) - M(t) (14) 

is the total available mass to be accreted by the NS, Mns(0 is the 
NS mass, and v e j re i(f) is the velocity of the ejecta relative to the 
NS 

In Eq. (TBI) . M e j(0) is the given initial mass of the ejecta (just 
at the beginning of the accretion process); we choose different 
values for it in Table|4] M(t) is the mass of the ejecta that is lost 
because it passes through the capture region of the NS. 

The actual mass accretion rate onto the NS, M accl -(f), is a frac- 
tion j] aca < 1 of Eq. ( fT3l , i.e. 



mass-radius relation of Belv edere et al.l (1201 2b . From the con- 
straint given by Eq. (TT~ST > we fix the binary separation a. We 
then proceed with the numerical integration of the accretion rate 
equations by requiring that Mns(0 = M cr it at t — Af accr , given by 
Eq. ( fT9l , from which we obtain the efficiency ?7 accl . 

It is interesting to analyze how the NS can accrete such a 
large mass, in some cases of the order of 47% of the early-SN 
material (see column 5 of Table @), since one could think that 
solid angles of ~ 50% between the early-SN material and the 
accreting NS are hard to obtain. Indeed, during the accretion 
process the NS is moving with high orbital velocities of the order 
of 10 8 cm s _I relative to the core progenitor (see column 7 of 
Table |4), and consequently travels effective arc-lengths several 
times larger than the circumference of the orbit (see column 8 of 
TableH. 

Assuming that the gain in gravitational energy of the accreted 
material into the NS can be released from the system leads to an 
upper limit of the luminosity 



M accr (f) = /7 accr M(f), 



(16) 



where 77 is the accretion efficiency onto the NS. So, there is an 
amount of material per unit time M oi „(f) = (1 - r) accl )M(t) not 
accreted by the NS. 

In Eq. (O, v O rb,Ns(0 = ^G(M pmg + M NS (t))/a is the orbital 
velocity relative to the SN core progenitor, a is the separation 
distance between the NS and the SN core progenitor, and 



VelW 



dr ei (t) r ej (f) 

— ; — = b 

dt t 



(17) 



is the expansion velocity of the early-SN material, where we 
have used r e j(?) = r em (f), given by Eq. dTTb - 

We have already mentioned that the power-law component in 
the spectrum of Episode 1 might be due to the accretion onto the 
NS companion. As this power-law component is present since 
the beginning of Episode 1 , we have fixed the value of fo, aC cr to 
be equal to the starting time of Episode 1 . This puts a constraint 
in the separation distance a of the binary, which under these con- 
ditions is given by 



a = r + fl cap (0), 



(18) 



where ro = r fJ (0) and R cap (0) are the radius of the early-SN 
ejecta and the capture radius of the NS companion at the begin- 
ning of Episode 1. In this case, ro ~ 1.75 x 10 9 cm, see Fig. [T2l 
The separation a is a function of the initial mass of the NS and 
of the SN core progenitor mass, as well as of the orbital velocity, 
through R cap . It is clear that the constraint given by Eq. ([T"8~t is a 
lower limit, since the accretion process onto the NS could have 
been triggered before by layers at lower densities (e.g. He). In 
such a case, the binary separation a could be higher. 

In addition to this constraint, we must take into account that 
the NS must reach its critical mass M cr j t at the beginning of 
Episode 2, since by that time the NS must collapse to a BH and 
emit the canonical GRB. This implies that 



Af a , 



611 
(1+z) 



349 s. 



(19) 



We show in Table|4]the parameters of the binary system lead- 
ing to IGC of the NS in a time interval equal to the duration of 
Episode 1 . We adopt an initial mass for the NS, M NS (0) = 1 .4M 
and, correspondingly, a NS radius of ^ns(O) = 12.3 km from the 



\E h (t)\ 



GM accl -(QM NS (f) 



(20) 



where we take into account the dependence of the NS radius with 
time, due to the increment of the NS mass by the accretion pro- 
cess. The self-consistent rad ius is computed at e ach time from 
the mass-radius relation of lBelvedere et al.l (120121) . 

The actual luminosity depends on the efficiency T] a i in 
converting gravitational energy into electromagnetic energy by 
some still unknown process. Since in our model we assume that 
the BB component of Episode 1 is due to the early-SN expan- 
sion, we estimate the efficiency r] m A from the assumption that 
\Eb\ is responsible for the power-law luminosity Lpl, namely 



77rad(0 



-PL 



\E„(t)\ 



(21) 



In Fig. \T5\ we show the evolution of the efficiency 77,3d in the 
first seconds of emission for the binary systems shown in Table 
[4] We assume a constant and isotropic power-law luminosity of 
Episode 1, L PL * 1.8 x 10 50 erg s _1 * 10~ 4 M s" 1 , as found 
from the spectral analysis. For all the cases, we obtain the same 
evolution of the efficiency with time, i.e. the curves overlap. 
This is due to the fact that we constrained all the systems to have 
the same initial NS mass and Af accr . 



8. Radio observations 



Zauderer & Be rgerl (1201 2l) report observations with the EVLA 
radio telescopes on several occasions between 11 and 16 July, 
at a frequency of 5.8 GHz. They found a radio source which 
brightened by about a factor of 1 .6, between 2. 1 and 7 days after 
the burst. The coincidence with the XRT position and the rising 
flux indicate that this is the radio afterglow of GRB 110709B. 
The position of the source is RA = 10:58:37.114, DEC = - 
23:27: 16 .760. We show in Fig . [16]the 5.8 GHz light curve pre- 
sented in Zaud erer et al.l (120121) . where there is evidence of a ra - 
dio bump. Following the work of lChevalier & Sod erberg ( 20101) . 
we have reproduced the plot of the peak spectral radio Luminos- 
ity per unit frequency as a function of time (days) at which the 
peak is produced for different SN associated with GRBs, includ- 
ing GRB 1 10709B (see Fig. [17]). We find that the radio emission 
of this source is higher than the ones associated with typical SN. 
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M prog /M 


M ej (O)/M 


p ej (0) (g cm- J ) 


?7accr 


AM accr /M ei (0) 


P (min) 


v'oi-b,Ns(0) (kms" 1 ) 


A^accr/ P 


a/R 


4 


2.7 


2.39 x 10 5 


0.92 


0.47 


0.52 


5.24 x 10 3 


11.14 


0.037 


5 


3.7 


3.27 x 10 5 


0.88 


0.34 


0.45 


5.84 x 10 3 


12.96 


0.036 


6 


4.7 


4.16 x 10 5 


0.88 


0.27 


0.39 


6.39 x 10 3 


14.71 


0.035 


7 


5.7 


5.04 x 10 5 


0.89 


0.22 


0.35 


6.91 x 10 3 


16.39 


0.034 


8 


6.7 


5.93 x 10 5 


0.91 


0.19 


0.32 


7.40 x 10 3 


18.00 


0.033 


9 


7.7 


6.81 x 10 5 


0.94 


0.16 


0.30 


7.87 x 10 3 


19.55 


0.032 


10 


8.7 


7.69 x 10 5 


0.96 


0.15 


0.27 


8.32 x 10 3 


21.04 


0.031 



Table 4. The massive star - neutron star binary progenitor of GRB 1 107 09B. M pl . oe is the mass of the massive star (in solar masses), M e j(0) is the 
mass of the ejected material in the early-SN phase (in solar masses), p c j(0) is the density of the ejecta at the beginning of the expansion, ^ accr is the 
efficiency of the accretion process onto the NS, AM accl . = M clil - M NS (0) is the total mass accreted by the NS before the collapse, P = 2^ra/v orb NS 
is the period of the binary, v ol -b.Ns(0) is the initial orbital velocity of the NS and At dca /P is the arc-length travelled by the NS during the accretion 
process in units of the length of the whole orbit and a/R Q is the binary separation (in units of solar radii). We suppose that the accretion process 
starts 5 s before the first trigger, i.e. fo. accr coincides with the time corresponding to the first datapoint in Fig. [5] 
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Fig. 15. Theoretical estimation of the efficiency ;/ rad given by Eq. l !2 1 1 
of the process to convert gravitational energy in radiation as a function 
of time. For this plot, we have assumed a constant and isotropic power- 
law luminosity of Episode 1,L PL « 1.8xl0 50 erg s -1 ra 10~ 4 M o s~'. We 
computed the values of the efficiency for the binary systems shown in 
Table|4] For all the cases, we obtain the same evolution of the efficiency 
with time, i.e. the curves overlap. The values of rj mA are always < 10%. 
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Fig. 16. X-Ray (black), radio (blue) and NIR (red, upper l imits) light 
curves of GRB 1 10709B. Taken from lZauderer et al.l dioilT) with kind 
permission. 



9. Conclusions 

GRB 110709B is a very peculiar source, since it is the first 
for which Swift BAT has triggered twice. Its Swift BAT light 
curve presents two well defined episodes, Episode 1 and Episode 
2. Episode 1 lasts 100 s and Episode 2 lasts 135 s. Particu- 
larly interesting is the fact that the X-ray observations started 
well before the second trigger. The light curve and spectrum 
of this source share similar characteristics with GRB 090618 
(Izzo, Rueda & RufBnill2012h. GRB 1010 23 dPenacchioni et alJ 
120121) and GRB 970828 dlzzo et alJfloil) . It has been recently 
shown that these GRBs which show such distinct emissions, 
Episodes 1 and 2, form a new family of GRBs described by the 
IGC paradigm dRueda & Ruffinil 120121: llzzo. Rueda & Ruffinil 
12012b . Within this scenario, the GRB originates in a binary 
system formed by a massive star on the verge of a SN and a 
NS close to its critical mass for the gravitational collapse to a 
BH. The compact core SN progenitor ejects material in the very 
early phases of the SN explosion, that is then accreted by the 
NS; this process is identified with Episode 1. The accretion 
process onto the NS brings it to the critical mass, leading to 
its gravitational collapse to a BH and emitting the GRB, iden- 
tified with Episode 2. Later on, we see a standard emission in 
X-rays, which we have called Episode 3. Several days after 
the burst, when it is present, we see an optical emission, asso- 
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Fig. 17. Plot of the peak spectral radio Luminosity per unit fre- 
quency versus the time at which the peak occurs, for different SN as- 
sociated to GRBs. The circles represent the SN emission associated to 
SN 2006aj (GRB 060218), SN 1998bw (GRB 980425) and SN 20031w 
(GRB 031203). The triangles represent the SN Ib/c for which there 
are radio observations, namely SN 2002ap, SN 1990B, SN 2008D, SN 
19941, SN 2009bb and SN 2003L. The squares represent the SN lib: 
SN 2008ax, SN 2001ig, SN 1993J, SN 2001gd and SN 2003bg. The 
red cross is the luminosity related to GRB 110709B afterglow. It is 
higher than the emissions of the other SN, considered "standard". 
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ciated to the S N (Episode 4). Following t he recent works on 
GRB 090618 dlzzo. Rueda & Ruffinil 12012) and GRB 970828 
dlzzo et al.l 1201 2|) . we here apply the IGC paradigm to GRB 
110709B. 

The redshift of GRB 110709B is unknown, so in Sec- 
tion |4] we use d four phenomeno logical methods to c onstrain 
it; i.e ., Grupe dGrupe et al.ll2007h . Amati (lAmatill2006h . Yone- 
toku dYonetokull2004h and the scaling of the X-Ray afterglow 
dlzzo et all l2012t iPenacchioni et alj 120121) . The first method 
gives an upper limit of z < 1.35. The second and third meth- 
ods give a lower limit of z > 0.6 and z > 0.7, respectively. The 
last method gives a precise value of z = 0.75, which lies within 
the range determined by three the above mentioned methods. We 
then fixed this last value as the redshift of GRB 1 10709B. 

The spectral analysis of Episode 1 is given in Section 13.11 
and [5] We find a value of the isotropic energy for Episode 1 of 
E^l = 1.42 x 10 53 erg (see Table [1). We fit the spectrum with a 
BB+PL model. The temperature of the BB component evolves 
with time following a broken power-law (see Fig. |5). The cor- 
responding radius of the BB emitter evolves in time following a 
power-law given in Eq. dTTT i and shown in Fig. [12] We asso- 
ciate this radius and the BB component to the evolution of the 
SN ejecta, while the power-law is associated to the accretion of 
the ejected material onto the NS companion. 

Episode 2 is analyzed in Sections 13.21 and [6] We find an 
isotropic energy of eF 1 = 2.43 x 10 52 erg (see Table [3). We 
interpret this episode as a canonical GRB and simulated its light 
curve and spectrum within the Fireshell model. We find at trans- 
parency a Lorentz factor T ~ 1.73 x 10 2 , laboratory radius of 
6.04 x 10 13 cm, P-GRB observed temperature kT P - CRB = 12.36 
keV, baryon load B - 5.7 x 10~ 3 , P-GRB energy of E P - GRB = 
3.44 x 10 50 erg, and a CBM mean density 76 p art cm" 3 . This 
value i s consistent with a "dark GRB", as cited in lZauderer et alJ 
d2012l) . The lack of detection of a SN emission for this particular 
GRB could be due to obscuration by the circumstellar dust in the 
host galaxy. 

The nature of the progenitor is discussed in Section[7] We in- 
dicate that it is a binary system formed by a massive evolved star 
on the verge of a SN explosion and a NS. We associate the ther- 
mal component of Episode 1 mainly with the early-SN evolution 
and the power-law component to the accretion process onto the 
NS. There is the possibility that also the accretion process has a 
thermal contribution. The energy due just to the thermal com- 
ponent is of the order of 10 50 erg, which is reasonable for the 
expansion of the early-SN ejecta. We perform all the necessary 
calculations to obtain the parameters of the binary system. For 
all our calculations we assume a fixed NS mass of 1.4M©. We 
compute the rate at which the early-SN material enters the cap- 
ture region, for given values of the SN core progenitor mass. 
From this material, only a fraction will be accreted by the NS, so 
we introduced an efficiency factor ?7 accl . As the power-law com- 
ponent is present since the beginning of Episode 1, we suppose 
that this episode starts at the same time forcer as the accretion 
process, namely, when the outermost shell of expanding ejecta 
reaches the capture radius R cdp of the NS (measured from the 
center of the NS). This puts a constraint on the separation dis- 
tance a of the binary. In addition, the NS must reach its critical 
mass and collapse to a BH at the beginning of Episode 2. This 
puts a constraint on the duration of the accretion process Af accl . 
By integrating the accretion rate equations with these boundary 
conditions we obtain the efficiency 77 accr . We summarize the re- 
sults in TableHJ for different values of the core-progenitor mass 
and the density of the early-SN ejecta. Assuming that the power- 



law radiation comes from the conversion of the binding energy 
of the accreted material onto the NS, we estimate the efficiency 
r] m d of this conversion process, which we show in Fig. [15]for an 
isotropic power-law luminosity Lpl ~ 1.8 x 10 50 erg s _I « 10~ 4 
M Q s _1 observed in Episode 1. For the parameters of the binary 
system shown in Table |4j we obtain values of ?7 ra d < 10%. The 
efficiency of the radiation mechanism can be even lower if some 
beaming or boosting is present. However, we did not address 
any such possible mechanism in this work. 

In section [8] we present the radio observations of GRB 
110709B with the EVLA radio te lescopes and the X-r ay, radio 
and NIR light curves taken from lZauderer et al. l (l20H . We no- 
tice the presence of a bump in the radio afterglow, at « 10 days 
after the burst. As GRB 110709B has been classified as an op- 
tically dark burst, we plotted the peak spectral radio luminosity 
per unit frequency as a function of time and compared it with the 
luminosities of typical SNe, to see if it was possible to find any 
coincidences that may indicate the presence of the SN in the ra- 
dio band. However, the luminosity we find is much higher than 
the ones of the standard SNe. 

We interpret, within the IGC paradigm, that GRB 110709B 
is a new member of the IGC family, in addition to GRB 090618, 
GRB 101023 and GRB 970828. 

A remarkable support of the above IGC paradigm comes 
from the observations of the X-ray afterglow emission of these 
systems. The X-Ray light curve is composed of an early steep 
decay, a plateau and a late decay. The analysis of the late decay 
of the afterglow luminosity has been identified with the cooling 
of the newly born NS, left by the SN explosion (Negreiros et al.l 
|20H . 
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